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ABS TRACT 


By analyzing CONREP data for oiler (AO) and ammunition (AE) ships 
replenishing destroyers (DD) and attack carriers (CVA) with regard to 
day and night operation and load requirement, it is shown through Chi- 
Square goodness-of-fit tests that alongside time data can be fitted to 
an Erlang or exponential distribution. In addition, several methods 
for estimating the shape, scale, and shift parameters for gamma and 
Erlang distributions are presented. 

Employing well known but analytically little used features of the 
replenishment operation, a computer simulation model is then formulated 
and programmed based on these distributions, and thus sensitive to 
changes in load requirements. Using the simulation as an experimental 
device, an example is run to demonstrate its use in estimating the time 
to complete the operation and conducting sensitivity analyses on the 


various load requirements. 
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I. INTRODUCTION 


The primary reason the United States Navy is capable of performing 
its regularly assigned mission as well as assist in global contingencies 
effecting United States' interests, is the ability to swiftly resupply 
combatant ships almost anywhere in the world. This is possible through 
replenishment at sea operations which render ships independent of 
costly, vulnerable land based supply facilities while effectively 
increasing combatant force strength, since ships are able to remain 
on or near their assigned stations during replenishment. 

Although very important, replenishment at sea is a hazardous 
operation. During replenishment, participant ships are extremely 
close increasing their possibility of collision and vulnerability to 
enemy attack. Even though vertical replenishment (VERTREP) using 
helicopters has helped to reduce the risk in some cases, all fuel 
transfers such as Navy Special Fuel Oil (NFSO), Aviation Gasoline 
(AVGAS), and JP-5, and most transfers of dry material such as ammuni- 
tion and stores, are necessarily done by the more dangerous connected 
replenishment (CONREP) method. Recognizing the importance of CONREP, 
the U.S. Navy is using a considerable amount of its resources to 
increase the efficiency and reduce the danger of this operation. 
Increased training of personnel, redesign or replenishment ships and 
equipment, and studies utilizing analytical and computer simulation 
models are efforts that have been and are being made in this endeavor. 

Operations Research studies of the replenishment at sea operation 
have been very few. However, McCoullough [4 | formulated two models of 


a simplified underway replenishment operation. One an analytical 
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infinite cycle model using multi-stage queueing theory and the other 

a non-cyclic computer simulation. In both models the distribution of 
service time at each replenishment ship was assumed to be exponential. 
Waggoner (7) and to some extent Patterson (6) also used the assumption 
of exponential service (i.e. alongside) times in replenishment at sea 
operations. Besecker C1 questioned the validity of this assumption 
and analyzing CONREP data, showed that the distribution of service times 
can be fitted with a gamma distribution. Categorizing the data as to 
participants and day and night operation, he contended that the data 
would also fit the more useful Erlang distribution in each case. This 
paper continues the investigation of replenishment at sea operations 


along the lines of reference (1) : 


12 


II. OBJECTIVES 


The most immediate gain in CONREP efficiency can be experienced 
through proper planning of the operation. However, studies of the 
operation have been largely ignored in the planning stage, since they 
lack the realism necessary to gain the confidence of operational 
commanders. Previous models (such as \4] and [7\) have not taken 
full advantage of the many known features of CONREP, for instance, 
arrivals at supply ships are not random but may be promulgated in 
advance by the Task Group Commander. Also, the quantity required by 
each combatant ship is known prior to the operation to some degree of 
accuracy rather than being a random variable. If the transit time 
between various supply ships is considered deterministic and known, 
the uncertainty in the system is a result only of the unknown alongside 
times, which should be dependent on the participants and the loads 
transferred. 

The first objective of this paper is to develop a more accurate 
estimation of alongside times for various combinations of CONREP 
participants by considering a given load requirement. Such conditional 
distributions would be applicable in analytical models, in computer 
simulations, and possibly in planning replenishment at sea operations, 
with more accurate predictions resulting in each case. 

The formulation and programming of a simulation model based on 
these alongside time distributions, and thus sensitive to changes in 
load requirements is the second objective. While incorporating the 


above known features, the simulation should be designed to enable its 
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user to easily adjust the arrangement of the ships involved, load 
requirements, and distribution parameters. In this regard, he could 
utilize it to test the suitability of various arrangements of combatant 


ships and the sensitivity of load changes. 
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III. SOURCE OF DATA 


Data used to obtain the distribution of alongside completion times 
were from replenishment at sea (RAS) reports 3180/1A and 3180/1B for 
the period December 1967 to April 1969. Commander Service Force, U.S. 
Pacific Fleet requires these reports be made to him immediately 
following underway replenishment by those replenishment ships under 
his cognizance. Each report includes the replenishment ship type 
and hull number, receiving ship type and hull number, total alongside 
time, quantity transferred in barrels or short tones as appropriate, 
UNREP speed, whether day or night, whether EASTPAC or WESTPAC, and 
whether or not there was simultaneous port and starboard UNREP in 
progress. Full instructions concerning the submission of these 
reports are available in COMSERVPAC INST 3180.3F of 6 November 1967. 

Alongside time is defined as the elapsed time from commencing 
approach from 300 yards astern to time of departure. Only data 
submitted by oiler (AO) and ammunition (AE) type ships involving 
destroyer (DD) and attack carrier (CVA) recipients were used, since 
a sufficient quantity was available only for these combinations. 
However, sufficient data should become available in subsequent years 


to make similar analyses possible for all combinations of participants. 


IV. METHOD OF DATA ANALYSIS 


In order to accomplish the objective of a more accurate estimation 
of alongside time distributions for various participants, data was 
first categorized as to type of ships involved, then separated with 
regard to day or night operation, and finally subdivided into sets 
based on quantity transferred. The distinction based on day or night 
conditions was considered necessary, since from personal experience 
and as shown by Besecker, SAN: operations consistently take longer 
than those conducted in daylight hours. The load intervals were 
selected on the basis of quantity of data available and uniformity 
of interval ranges. The latter was desirable in order to establish 
load intervals necessary for the computer simulation. 

Data in each category was analyzed by constructing histograms of 
the relative frequency with which the alongside replenishment time 
fell within ten minute intervals. The last interval is always a half 
open interval from its selected lower bound to infinity. In each case 
the histogram suggested that the distribution of the completion times 
for the respective load intervals could be fitted to a gamma or 
exponential distribution presented in Appendix A. This was expected 
since tie distribution of alongside times over all intervals had been 
fitted to a gamma distribution in fel: 

Figure (1) is the histogram constructed for the case in which an 
oiler (AO) was transferring 1000-1099 barrels of fuel to a destroyer 


(DD) during daylight hours. This is presented as a typical example. 
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It was appropriate to apply a Chi-Square goodness-of-fit test to 
the null hypothesis 


He ea enh cae) 


against all alternatives, where P, is the probability of finding a 


i 
completion in the fea interval with respect to the distribution in 
question, P,;,, is the relative frequency of the data falling within the 
same interval, and k is the number of intervals. The probability P 
was calculated using Besecker's computer subroutine GAMDIS with the 


gamma parameters estimated from the sample mean and variance using 


the relationships 


7 
Aum ce 
Yr 52 
\_e 
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= 2 
where T is the sample mean and s' the sample variance. 


The Chi-Square test consisted of calculating the statistic 


MP io 


K 
oe \ (X;-nP, ve 
i=1 
where is the number of datum which fell within the ia ten minute 
interval and n is the sample size. In accordance with statistical 
theory 3 , if H is true, then the random variable Q has an approxi- 


mate Chi-Square distribution. The hypothesis was accepted if 


2 
O<% C95, tad). 


The degrees of freedom was necessarily k-3 since two parameters were 


estimated from the data. Many authors caution that nP,, ee 2 5 dle 
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should be greater than five, since Q has only an approximate Chi-Square 

distribution. This warning was strickly adhered to, which necessitated 

combining two or more time intervals when the quantity of data in a load 
interval was small, 

In the cases involving attack carriers, the hypothesis was accepted 
for all load intervals. Since the Erlang distribution is more ameanable 
to analytical and computer simulation models and only differs from the 
gamma distribution in that tr is restricted to be a positive integer, each 
was rounded to the nearest positive integer, the parameter *) adnmmeel 


according to the expression 


, 


Hi | > 


and a similar test conducted. Again, the hypothesis was accepted in 
each case. Results of these latter tests are given in Tables I, II, 
III, and IV in Appendix B. Figure (2) shows a typical example with the 
probability density function superimposed on the relative frequency 
histogram for a case involving AO-CVA participants. The figure illus- 
trates the method of combining intervals consistent with the greater 
than five rule given above. 

When the above method was used for cases involving destroyer type 
ships, the hypothesis was rejected in each case, Re-examination of the 
histograms for these cases suggested that a shift of the time axis by 
an amount a would lend to an acceptance of the gamma distribution. 
Several approaches were taken to estimate the distribution parameters 
©. , and a. Although many were not successful in this study, all are 
presented for reference in future data analyses. 


Noting that the first three centeral moments of the gamma distribu- 


tion given by 
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could be used to estimate the parameters, the following relations were 


derived: 
A 6 
oe (2) 
ae. 
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N gt 
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where T is the sample mean given by 
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2 oe ; 
s’ is the sample variance given by 
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and a is the third centeral moment given by 
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By rounding the shape parameter r to the nearest positive integer and 
adjusting the scale and shift parameters, the following expressions 


were obtained: 
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where (x)* is the nearest integer to x. Using this method in several 
cases, = was found to be unity in each case suggesting an exponential 
distribution. Figure (3) shows a typical case with the exponential 
distribution superimposed on the histogram. The hypothesis was rejected 
for these cases as clearly suggested by the figure. 

Concern over the poor results and the fact that observed times were 
less than a, prompted another investigation into the method of obtaining 
the shift parameter. Hypothesizing an exponential distribution and | 


constructing the likelihood function 


LC Age) = Nema 2 for ai min T 
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the maximum likelihood estimator for a was obtained as 
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from Eq. (1). However, for all cases examined using this method, the 
hypothesis was rejected on the basis of the Chi-Square goodness-of-fit 
test as illustrated in Figure (4). 

Concluding that rounding to obtain an Erlang distribution was 
inappropriate when os ré 1, an attempt was made to use the maximum 
likelihood estimate for a under the assumption of the gamma distribution. 
The likelihood function in this case is 


n 
Tn a) (T,; - a) 
i=l 


n 
r-1 

L(a) = d W(T, - a) = for atmin T 

U= \ i ASUS mM i 

0 otherwise, 


resulting in the maximum likelihood estimate of a as 


again. Then referring to the expressions obtained using the method of 
moments and adjusting the other parameters to a, the following equations 
were derived: 
~ 2 
7 Clea) 
S 
A = 
d= as 


2 
S 


Gamma distributions suggested by this third method did not fit the data 
as shown in Figure (5). 

In another attempt to use the method of moments to estimate the 
distribution parameters, no adjustment was made to the estimates 
given by Eqs. (2), (3), and (4). But this also proved unsatisfactory 
as illustrated in Figure (6). 

Convinced that the gamma distribution was appropriate from the 


shape of the histograms based on ten minute intervals from zero, 
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extensive probability tables were generated using Besecker's subordinate 
GAMDIS. Then by scanning these tables for similar patterns of the 
relative frequency ratios of time intervals compared to the probability 
that a completion time would be within the same interval, an Erlang 
distribution was selected that appeared to fit the data. A Chi-Square 


goodness-of-fit test was conducted and the hypothesis H accepted if 
Z 
Q X (.95,k-1). 


This method proved quite fruitful, resulting in the acceptance of 
the hypothesis associated with an Erlang distribution in each case. A 
typical case is shown in Figure (7). Results of these tests are given 
in Tables V, VI, VII, and VIII in Appendix B with T calculated using 


the expression 


ee} > 
Il 

> >|n > 
+ 
fo > 


The estimate a was chosen as the lower bound of the first interval in 
which data was observed. Figures (12) through (56) in Appendix C 
illustrate the fit of the probability density functions and the histo- 


grams for each case involving destroyer type ships. 


28 


OS1 OOT US 0 
YO - 


PeAIesqO seM e7ep YoOTyM ut 





O°O€ = [BAT@IUT 3s2TZ eYI Jo punoqg temo, = e 

| OC ° 

Z¥ ; | 

OOT a a7 =e | | 

7 ! f 

a is 

= a = Vv — , 
y Ue" 


"eased qd ‘SA OW [eoOTdA e@ X07 weasgojsty AouUenbailzZ SATIBTAI sAzTQOedSOI 
94 uo pasodut poyjew gutuueds Butsn uotzounjz AQtsuep Ajt{tqeqoid ewues pajses3ng ‘/{ sAan3iy 


OfjJe1 ADSuonbs1y 
2° 


V. THE COMPUTER SIMULATION MODEL 


As stated in Naylor [5] the purpose of simulating a large operation 
such as replenishment at sea is to conduct situational experiments that 
would ordinarily be too expensive and/or cumbersome to perform physi- 
cally. Because of the many conceivable combinations of participant 
ships and variations in quantities transferred between them that one 
might want to examine, physical experimentation would be unreasonable 
here. Therefore, based on the distributions of the previous section, 
a Monte Carlo type digital computer simulation of the replenishment at 
sea operation named SIMA was devised and programmed in accordance with 
the second objective. 

SIMA was designed to accommodate a system of from one to four 
replenishment ships, composing an Underway Replenishment Group (URG), 
with each ship having one or two replenishment stations on the port or 
port and starboard sides. A maximum of twenty combatant ships can be 
in the queue at each station. Furthermore, six types of combatant 
ships can be in the system. This capacity is well beyond what one 
would normally expect during an actual replenishment operation, 
therefore overload of the model is not expected. 

Drawing from personal experience, the model was formulated as 
realistically as possible, thereby taking advantage of the known 
factors of a replenishment at sea system. All combatant ships are 
initially distributed among the various stations of the supply ships, 
with each individually characterized by the load it requires from each 
replenishment ship. Combatant ships proceed sequentially through the 


System and upon arrival at a replenishment ship it enters the smaller 
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(if more than one exists) queue. Each ship receives the services of 
every supply ship and, having done so, is ejected from the system. 
Although the combatant ships do pass through the system in a specified 
manner from their initial positions, these conditions closely resemble 
the ability of a Task Group Commander to formulate the operation prior 
to its commencement, utilizing his knowledge of the approximate loads 
transferred among the participants. 

Using the simulation as an experimental device, one is able to 
specify the participants, number of stations, initial sizes of the 
station queues, and quantities transferred between the respective 
ships, and then periodically observe the state of the system. These 
observations are made at variable time increments which are the random 
alongside completion times generated as per Appendix D. An observation 
is made and the system advances whenever any one of the URG ships has 
completed a replenishment. Each observation consists of the time since 
the start of the operation and the state of the system. The state of 
the system at these times is the number of ships, including the one 
being served, at each replenishment station and is represented by an 
n dimensional vector, where n is the number of stations in the system. 
The ‘i component of the state vector being the number of ships at the 
co station. 

This model was programmed in Fortran IV and tested on the IBM 
model 360-67 digital computer at the U.S. Naval Postgraduate School, 
Monterey, California. A printout of the program is presented in 
Appendix E. Because of the digital nature of the simulation, ship 


types were given the following identification numbers: 
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Oiler - l 
Ammunition ship - 2 
Destroyer - l 
Attack carrier - 2 


Since all the parameters derived in the previous section are given in 
data statements, they may be revised by changing the statements. In 
using this technique, a user need only supply the following information 


to conduct replenishment at sea experiments: 


1. Type and sequence of supply ships 

Day or night operation 

3. Type, initial position, and load requirement of the combatant 
ships. 


ho 


This information is provided by a data deck. The first m(m - 1,2,3, or 


4) cards represent the m supply ships, each card having the following 


format: 
Columns 1-3 Type of replenishment ship, i.e. AO, AE 
Columns 4-5 Number of stations available. (1 or 2) 
Columns 6-7 Supply ship identification number 
Columns 8-10 Blank 
Column 11 Must be either 0,1, or blank. 


0 or blank implies that there is at least 
one other replenishment ship in the system and 
another card must be read. 

1 implies no more supply ships are in the 
system and the program continues to the next 
statement. 

Columns 12-80 Blank 


The order in which these cards appear designates the sequence in which 
combatant ships pass through the system. The m+1°° card specifies day 
or night operation and has the following format: 
Column 1 Blank 
Column 2 Must be either 1 or 2, 
1 designates day 
2 designates night 
Columns 3-80 Blank 


The remaining cards give information concerning the combatant ships, 


with the cards arranged as these ships initially appear in the system. 
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The remaining cards give information concerning the combatant ships, 

with the cards arranged as these ships initially appear in the system. 
The initial condition of each station's queue for each supply ship is 
given in accordance with the order of the supply ship cards, with the 
port station designated first. Each card represents a combatant ship 


with information given in the following format: 


Columns 1-3 Ship type, i.e. DD, CVA 
Columns 4-5 Ship identification number. 
Columns 6-11 Load required trom first supply ship. 


Columns 12-17 Load required from second supply ship. 
Columns 18-23 Load required from third supply ship. 
Columns 24-29 Load required from fourth supply ship. 
Column 30 Must be either 0, 1, or blank. 

O or blank indicates that there is 
at least one other combatant ship in 
the queue for the station in question 
and another card must be read, 

1 indicates that no more ships are 
initially in the queue in question. 

Columns 31-80 Blank 


Figure (8) shows a hypothetical initial arrangement of the supply 
and combatant ships, illustrating the possible routes of the combatant 
ships. As an example, this system was simulated using SIMA in day and 
night operation. Figure (9) shows the data cards submitted to set up the 
System for the daylight case. Since only two supply ships are in the 
example, the results may be represented graphically as in Figures (10) 
and (11). The origin represents the starting time of the operation. 
Each dot represents a combatant ship having completed replenishment at 
either the AO or AE. The length of the line segment between two 
consecutive dots indicates the elapsed time in minutes between two 
consecutive completions of the replenishment operation. For example, 
in Figure (10) a ship departed the AO at approximately 60 minutes, a 


Ship left the AE at approximately 80 minutes, a ship left the AO at 
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Hypothetical initial arrangement of supply and 
combatant ships. 


Figure 8. 





Supply ships as designated 


Destroyer type ship (DD) 


\ Attack Carrier (CVA) 
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Figure 10. Graphical representation of simulation 
results for the example system during 
daylight operation. 
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Figure 11. Graphical representation of simulation results for the 
example system during night operation. 


a7 


approximately 100 minutes, a ship left the AE at approximately 100 
minutes, and so on, regarding the commencement of the replenishment 
operation O hours and 0 minutes. Knowing the initial arrangement of 
the ships, one is able to reconstruct the movement of each combatant 
ship through the system and its departure time from the supply ships. 

Upon running the simulation fifty times for the example, it was 
found that the expected (average) time to complete the replenishment 
operation during daylight hours was 598.62 minutes while night operations 
required 635.38 minutes. Figures (12) and (13) are the histograms for 
the frequency ratios of total replenishment time of all combatants for 
these experiments. 

We see from the data cards that the destroyer, initially the second 
ship in the queue at the port side of the AO, requires 648 barrels of 
fuel oil and 2200 short tons of ammunition. When this was changed to 
1100 barrels and 3000 short tons and the simulation rerun, the expected 
time to complete the operation was 605.54 for day and 638.10 for night 
conditions. Sensitivity analyses along these lines could easily be 
accommodated by the simulation model. 

SIMA = hie “slo be useful in planning replenishment at sea 
operations. By tracing the path of particular ships in simulation, 
one could detect possible bottlenecks in a proposed replenishment 
system. Then by repositioning one or more ships in the formation and 
conducting additional experiments, one could plan a more efficient 


arrangement to be used in actual replenishment operations. 


38 


Frequency ratio 


Frequency ratio 


40 


. 30 


. 20 


.10 


.00 


Figure 12. Rela 


tive frequency histogram for total 


replenishment time of all combatants from 
fifty simulation runs for day operation. 
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Figure 13. Relative frequency histogram for total 
replenishment time of all combatants from 
fifty simulation runs for night operation. 
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VI. SUMMARY 


In this paper we have shown that an Erlang or exponential distribu- 
tion may be fitted to the CONREP alongside time data involving AO, AE, 
DD, and CVA type ships for each of several load intervals. It is 
possible for these distributions to be used in analytical models 
einnploying Laplace transforms that would be sensitive to changes in 
load requirements of the combatant ships. Furthermore, several methods 
for obtaining the parameters for the gamma and Erlang distributions are 
presented for reference in future data analyses. 

The second part of the study utilizes the distributions in a compu- 
ter simulation model. Using the simulation as an experimental device, 
one can easily and realistically investigate the suitability of various 
CONREP situations confident that the results will reflect one's know- 
ledge of the approximate loads required by each of the combatant ships. 
It was demonstrated that SIMA could be used to estimate the time to 
complete the operation and investigate the sensitivity of certain load 
requirements, and suggested that the model be used in planning underway 


replenishment formations, 
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APPENDIX A 


The Gamma, Erlang, and Exponential Distributions 


The gamma probability density function is given by 


£ (x) a eo a x 
(a) for x20 
0 otherwise, 


where r is the shape parameter and ) the scale parameter. In shifting 
the function by an amount a the random variable x would be replaced by 
x-a in the above expression. 

The Erlang and exponential distributions are special cases of the 
gamma distribution. If r is a positive integer, the probability 
density function may be written as 


rn-l - x 


Coe A x e 
(n-1)! for x2 0 


0 otherwise, 


and is that of the Erlang distribution. If r is unity, the distribu- 


tion is the well known exponential distribution. 
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APPENDIX B 


Tables 


Table I. AO vs. CVA - day 


load int. T + r Q X95. k-3) 
G-3999 85.0 7096 8 ZeOLS ae 
4000-4999 Was aCe 24 2576 3.84 
5000-5999 yee le? 14 3.444 a. o9 
6000-6999 2 eZ .074 2 1,497 Oot 
7000-7999 Wee .050 E ‘eGo ooo 
8000-8999 39.5 . 156 22 2a 5264 
9000-9999 126.3 047 6 i.2103 a, 99 
10000- Lass .060 4 2.584 eee 
Table II. AO vs. CVA - night 
load int. T ) r Q 4.95 k-3) 
0-5999 one .068 8 1.007 Ieoe 
6000-7999 Poe 088 ily 1.154 eo 
8000-9999 P39 .063 9 0.856 oes 
10000-11999 141.8 . 136 13 0,732 Da 
12000- oA 2 .056 11 02635 Deo7 
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load int. 


0-9999 

10000-19999 
20000-29999 
30000-39999 


40000- 


load int. 


Oe oS 

20000-29999 
30000-39999 
40000-49999 


50000- 


Table III. 
- A 
c A 
74.7 rOo2 
105.2 114 
Looe? Pere 
Looms 141 
229.4 07? 
Table IV. AE vs. 
a A 
T D 
oe 084 
16225 20 
168.8 .076 
21627 ls 
7 oral 094 


AE vs. 
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eZ 


18 


ye 


18 


CVA 


BS 


24 


24 


CVA - day 


2.486 
2.624 
Sb, 
Zona0 


0.067 


- night 


2S 
0,541 
1.402 
0.649 


05623 


10.95, k-3) 
3.84 
3.84 
5.99 
3.84 


3.84 


2 
“X(.95,k=3) 
3.84 
3.84 
3.84 
3.84 


Sago 


Table V. AO vs. DD - day 


A A A A 
load int. t T D 15 a Q 0.95 ,k-1) 
0-399 41.3 40.0 . LOO 2 20 4.558 ba LO 
400-499 44.9 B10) AO) . LOO 5 20 DOD 9,49 
500-599 49.8 50.0 .098 2 30 ie o73 9.49 
600-699 43.1 40.0 . LOO Z 20 1.104 9.49 
700-799 At 51.4 . 140 3 30 Bey 2 eo | 
800-899 ee) Do . 060 2 20 eo te 12760 
900-999 32.20 bo, 0 E20 3 30 Tres 9.49 
1000-1099 49.6 Do). 0) 2 3 30 4.413 14.10 
1100-1199 5340 50.20 . 100 3 20 Dao 2o 7 o.0 
1200-1299 53.4 50.0 . LOO 5 20 Sr 167 P2760 
1300-1399 56.4 Goo . 140 5 30 e339 I eou 
1400-1499 54.1 550) 040 1 30 kod PZancU 
1500-1599 6226 58D 070 Z 30 4.637 12.60 
1600-1699 Go.) 6o5) . 140 5 30 4,189 12.60 
1700-1799 6.3 60.0 . LOO S 30 3.306 9.94 
1800-1899 64.1 60.0 . LOO 3 30 P2509 Pio 
1900-1999 G/23 70.0 . LOO 3 40 2.439 9794 
2000-2099 67.4 40.0 . LOO 3 40 57073 Ele 0 
2100-2199 O5e2 650 040 l 40 2.306 9.49 


2200- 66.3 70,0 . LOO 3 40 4.925 Pe 
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Table VI. AO vs. DD - night 


a A n h A 
load int. T T ee r a Q "(4.95 K-21) 
0-599 65.6 6725 .080 3 30 5.626 Pi 10 
600-699 6340 67.5 .080 3 30 1.50 9.49 
700-799 LOZ 6745 .080 3 30 Owm2i53 7 von 
800-899 78.6 Gia . 080 3 30 1.784 7.81 
900-999 67.8 67.5 .080 3 30 0.318 7 oe! 
1000-1099 61.3 67.5 . 080 3 30 1 «a3 9.49 
1100-1199 G2, 1 67.5 . 080 3 30 2912 7 3On 
1200-1299 56.1 Ga 0 . 080 2 40 Te Dino 
1300-1399 67.6 65.0 080 2 40 8.014 9.49 
1400-1499 71.6 67.5 . 080 3 30 los rac 
1500-1599 Do. 2 65.0 . 160 4 40 ao jzoe 
1600-1699 74.4 155 . 080 3 40 0.549 7 388 
1700-1799 7528 Jo . 080 3 40 1.854 7.61 
1800- 76.7 7 5¢5 . 080 5 40 4.563 9.49 
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Table VII. AE vs. DD - day 
load int. 7 t > r a Q (095 k-1) 
0-999 42.8 46.0 Zo 9 10 3,504 oa 
1000-1999 63.3 60.0 . LOO 4 20 1 [0 9.49 
2000-2999 7720 Sone 080 4 30 3. oe LOR 
3000-3999 STG 106.6 . 060 4 40 19 14.10 
4000-4999 110.5 105.4 044 Z 60 1.524 iS 
5000- 127.5 120.0 .080 4 70 5 0G 12.60 
Table VIII. AE vs. DD - night 
= A A A A 
load int. iL T r a Q ‘\rosmiee 
0-999 45.6 40.0 . 100 3 10 5.634 Jaa 
1000-1999 Te 67.5 080 3 30 io) oa 
2000-2999 70:9 S755 080 3 50 Ze lS oe 
3000-3999 104.7 — 1007.0 .040 Z 50 3. 600 O39 
4000- ie. 7 eo .080 3 70 eee 12,60 
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APPENDIX C 


Figures 


This appendix contains figures illustrating the fit of the 
probability density functions to the histograms for all cases 


involving destroyer type ships. 
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APPENDIX D 


Generating Random Observations from an Erland Distribution 


The most important operation in a Monte Carlo type computer 
simulation is the generation of appropriate random variates. In this 
study, it was necessary to obtain random variates from the Erland 
distribution. Methematically the Erland distribution is a convolution 
of r independent identical exponential distributions, i.e., the 
distribution of the sum of r exponentially distributed, independent 


random variables whose density is given by 


f(y) =e, yo. 


Thus, we are able to generate Erlang random variables by summing r 


independent exponentially distributed variables Yy> Yo» oS siens, 


~ 


each with expected value 1/ . Doing so we obtain that 


r 
x -\ Y; 


i=] 


has the Erland distribution with parameters r and , i.e. with density 


r r-l -Ax 
x 


aX e 
£ = A ee ee 
a (eel): ! for = 20 
0 otherwise. 


The question now is how to obtain exponentially distributed random 
variates? The cumulative distribution function for the exponential 


distribution is given by 


a3 


Cea neon ay 0. 


It is well known that F(Y), and by symmetry 1-F(Y) also, is a uniform 
random variable over the interval (0,1) regardless of what the distri- 


bution of the continuous random variable Y is. This is equivalent to 


setting 


where W is a random variable from the uniform (0,1) distribution. 


Therefore, taking the natural logarithm of both sides 
-AY = ln(W) 
which yields 


1 
> n(W) 


Thus, if W is a uniform random variable over (0,1), Y will have the 


desired exponential distribution. Substituting into Eq. (1) we get 


1 rE 
X =- —In(Tr W,). 
i 
ell 
There are several methods available for generating uniform (0,1) 


random variables. The IBM Scientific Subrouting RANDU was used in 


the study. 
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APPENDIX E 


Computer Program 


This appendix contains a printout cf the computer program SIMA. 
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GENTRATE ALONGSIDE TIMES FOR ALL INITIAL CUSTOMERS 
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GENERATE ALONGSIDE TIMES FOR THE SHIPS INITIALLY BEING SERVED. 


CALL GAMMACIX »NAAs AAAy XBy IX) 


OOW 


GO TO 21 


he Toe 
NNISEQe 2e ANDeo lL Neo EQe 1) 


es 
NUE 


: 


XQ 
IF 
27 CONT 


XB=XB4+LTAC(L 


VECTOR AND TIME. 


PRINT OUT INITIAL CONDITION 


DO 30 JI=1,!I 


UUW 


28 


e2eANDeMJeEQe1) GO TO 


NIQ(JI,MJ) 
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